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ABSTRACT.  Multicommodi ty  formulations  of  the  travelling 
salesman  problem(TSP)  use  usually  node  commodities  that  have 
to  be  send  from  a  given  source  to  the  various  nodes  of  the 
problem. An  alternative  multicommodity  formulation  for 
the  TSP  is  to  consider  the  arcs  of  a  solution  to  the  TSP 
as  commodities. 

1.  INTRODUCTION  Many  f ormultations  for  the  travelling 
salesman  problem  have  been  developped. These  include 
n-commodity  formulations  ( Wong[ 6] , Claus[ 1 ] )  and  2-commodity 
formulations  ( Finke, Claus, Gunn[ 3 ] , Finke[ 4] ) . 

The  formulation  of  the  TSP  described  in  this  paper 
uses  first  n2  commodities  to  describe  the  problem  and 
then  n6  to  proof  certain  theorems  about  the  solutions 
satisfying  the  required  constraints.lt  will  for  example 
be  shown  that  a  convex  combination  of  two  assignment  problems 
will  only  be  feasible  if  the  solution  is  also  a  convex 
combination  of  hamiltonian  paths. 

2^  THE  PROBLEM  AND  FORMULATION.   The  problem  we  consider 
is  one  of  finding  a  minimum  cost  Hamiltonian  path  from  one 
vertex  s  to  another  t  in  a  given  weighted  directed  graph. 

The  programming  approach  involves  introducing  a  variable 
Xij  for  each  pair  of  nodes,  that  takes  on  value  if  the  edge  i 
to  j  is  on  the  path  corresponding  to  the  given  choice  of 
variables  and  0  if  it  is  not.  (In  the  corresponding  polytope 
these  variables  are  constrained  to  be  between  zero  and  one.) 

The  continuity  constraints,  for  all  i, 

SUM(i)   Xij  =  1,  SUM(j)   Xij  =1, 

SUM(j)   Xsj  =  1,  SUM(j)   Xjt  =1, 

SUM(i)   Xis  =  0,  SUM(j)   Xtj  =0, 

insure  that,  for  an  integral  solution,  exacly  one  arc  enters 
and  leaves  each  vertex. 
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There  are  several  ways  to  add  additional  constraints  in 
order  to  guarantee  that  an  integral  solution  corresponds  to  a 
Hamiltonian  tour. 

The  Dantzig,  Fulkerson  and  Johnson[2]  approach  was  to 
require  that  for  every  nontrivial  subset  V  of  the  nodes  ( 
other  than  s),  we  have 

Xij  >=1 
for  all  i  not  in  V  and  j  in  V.  This  insures  that  one  can  get 
from  the  source  to  any  set  of  nodes  in  a  tour  corresponding 
to  an  integral  solution. 

In  {Claus[l])  n-commodities  one  corresponding  to  each  node 
were  used  to  define  the  subtour  elimination  polytope,and 
a  more  compact  two-commodity  versions  of  the  TSP  can  be 
found  in  ( Finke, Claus, Gunn[ 3 ] , Finke[ 4] ) . 

We  now  introduce  a  different  approach.  The  commodities 
kp 
X   now  to  be  send  from  the  source  are  the  arcs  (Xkp)  . 

kp 
Thus  for  each  arc  Xkp  we  define  variables  Xij      (which 

kp 
correspond  to  the  flow  of  commodity  X   on  the  arc  i,j)  with 
the  following  properties: 

kp  kp 

1.1  Continuity  at  i:  SUM(j)  X j i    =SUM(j)   Xij 

for  i  other  than  s  or  t; 

kp  kp 

1.2  Flow  from  s  to  k:  SUM(j)  Xsj   =Xkp,  SUM(i)    Xik   =Xkp; 


kp  kp 

1.3  Flow  from  k  to  t:   SUM(j)  Spj   =Xkp,  SUM(i)    Sit  =Xkp; 

kp       kp 

1.4  Continuity  at  i:    SUM(j)   {  Sij   -    Sji  J=0; 

kp       kp    kp 

1.5  Idendity       :        Aij   =   Xij  +  Sij; 

kp 

1.6  Capacity        :        0  <=   Aij    <=  Xij, 

kl  Ik 

SUM(l)   Aij  =Xij,  SUM(l)   Aij  =Xij; 


qr     kp 
1.7  Symmetry        :         Akp  =  Aqr; 
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kpqr 
A  second  level  of  commodities( variables)   A 
is  introduced  with  similar  types  of  constraint  as  the  above 
exept  that  the  quantities  leaving  the  source  and  transiting 

kp 
each  node  are  now  equal  to  Aqr      instead  of  Xkp,and 

kpqr  kpqr    kp 

2.1  Continuity  at  i:     SUM(j)   A j i    =SUM(j)   Aij   =  Aqr 

for  i  other  than  s  or  t; 


kpqr     kp  kpqr    kp 

2.2  Flow  from  s  to  k:  SUM(i)  Asi     =Aqr,SUM(i)  Aik  =Aqr; 

kpqr     kp 

2.3  Capacity:  0  <=   Aij    <=  Aij, 

kpqr      qr 
0  <=   Aij    <=  Aij  , 

kpqr   qrkp 

2.4  Symmetry:  Aij   =  Aij; 

kpqr    ijqr    ijkp 
Aij   =  Akp    =  Aqr 

Before  considering  the  flow  behaviour  for  a  convex 
combination  of  two  assignment  solutions  let  us  define  the 
following  terms: 

A  'Perfect  Exchange  Node' (PEN)  is  a  node  where  all 
flows  comming  into  a  node  on  one  arc  must  leave  the  node 
together  on  another  arc. 

A  'Regular  Flow'(RF)  is  a  flow  from  the  source  such  that 
along  an  RF  path  all  arc-commodities  progress  together 
along  such  a  path  Fig  1). 

A  'Regular  Flow  Path'(RFP)  is  a  path  with  RF. 

si          si  si          si 

Asl  >0  A12  >0  A23  >0      A34  >0 

12          12  12          12 

Asl  >0  A12  >0  A23  >  0     A34  >0 

(s) >(1) >(2) >(3) >(4) 

23           23  23          23 

Asl  >0       A12  >0  A23  >0      A34  >0 

34           34  34          34 

Asl  >0       A12  >0  A23  >0      A34  >0 

Fig  1 
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3. TWO  ARC  THEOREMS  MID  PROOFS 


Theorem 


Proof 


l:If  a  solution  has  an  RFP  of  length  h<n  then 
we  cannot  have  two  arcs  back  from  h  towards 
the  RFP. 
: Consider  the  RFP   shown  below  (Fig  2)  and 
assume  the  contrary  i.e  node  h  has  two  arcs 
back  towards  the  RFP.  In  this  case  the 

kpqr 
commodity  A      (were  hp  and  hr  are  the  two 
arcs  leaving  node  h)  cannot  leave  node  h, 
violating  constraints  2.1. 


(s) 


>(k) 


Fig  2 

Theorem  2: If  we  have  an  RFP  of  length  m  to  node  h 

and  node  h  is  not  the  sink  then  we  have  an  RFP 
of  length  m+1. 

Proof     : Consider  the  RFP  in  Fig  3 . Now  assume 

that  one  of  the  arcs  leaving  node  h  (h,f) 
has  a  flow  value  of  zero  for  commodity  (k,p) 

kp 
Ahf=0  .This  requires  that  this  commodity  as  well 
as  the  other  commodities  on  the  RFP  are  positive 

ed 
on  the  arc  (h,g)  (  Ahg   >0,  ed   RFP ). Due  to 


hg 
constraints  1.7  we  also  have  Aed  >0,  hg 
thus  prooving  the  theorem. 


RFP, 


(s) 


>(h 


Fig  3 


Theorem  3: If  we  have  a  regular  flow  of  length  n-1 
then  we  have  a  regular  flow  of  length  n, 
where  n  is  the  number  of  nodes  of  the  graph. 

Proof     :Consider  an  RFP  of  length  n-1  (Fig  4).  Node  n 
must  have  two  arcs  to  this  path. But  this  would 
require  that  the  path  is  not  an  RFP. 


(s). . .(k)--->(pj. . .(q)--->(r). . .(n-1 )--->( t) 


kp 
Anr  =Xkp 


Fig  4 


Theorem  4: A  path  from  s  to  t  cannot  be  an  RFP  unless 
it  is  part  of  an  assignment  solution. 

Proof     : Assume  we  have  an  RF  path  from  s  to  t  of  length 
h<n.Due  to  the  flow   constraints  (1.2)  we  must 
have  at  least  one  arc  connected  to  the  other 
nodes  not  on  the  RFP (Fig  5).  Assume  that 
this  node  is  node  r.Then  the  other  inward  arc 
into  node  r  must  have  positive  values  for  all 
the  commodities  on  the  RFP.  Thus  we  have  an 
RFP  path  going  through  arc  (q,r).This  path 
cannot  go  back  to  any  nodes  on  the  original  RFP 
path  and  is  thus  a  RF  cycle. If  there  are  more 
nodes,  then  one  of  those  must  have  an  inward  arc 
from  the  RFP  or  the  RF  cycle  etc 


(s)--->(l)---(2)--->(3)--->(4)--->(5)--->(6)--->(t) 

(f) >(io) 


34 
A912  =X34 
3  4kp 
thus     A912  >0 


( 


2)<---(ll) 


Fig  5 


-6- 

The  previous  theorem  showed   that  RFP ' s  are  only 
possibe  on  assignment  solutions . This  implies  that 
a  feasible  solution  to  the  arc  formulation  model 
is  a  convex  combination  of  RFP ' s . 

Theorem  5: A  convex  combination  of  two  RFP ' s  must  also  be  a 
convex  combination  of  two  Hamiltonian 
Paths, otherwise  constraints  (1.2)  cannot  be 
satisfied. 

Proof     : Assume  the  contrary  i . e . the  two  RFP ' s  consist 
of  zero  or  only  one  Hamiltonian  path. Since 
the  RFP'  carry  only  arcs  from  their  own  solution 
an  RFP  can  only  be  a  hamiltonian  path  in  order 
to  satisfy  constraints  (1.2). 


(s)--->(l)---(2)--->(3)--->(4)--->(5)--->(6)--->(t) 


10,7 
A   =0 


10,7 
A   =0 


>(9) >(\0)<- 


>(12)<- 


(11) 


10,7 
A  =0 


10,7 
■A   =0 


Fig  6 


Theorem  6: If  n-1  is  even  and  the  solution  for  n-1  nodes 
has  to  be  a  convex  combination  of  hamiltonian 
paths  then  it  is  also  true  for  n  nodes. 

Proof     : Since  n  is  odd  we  must  have  at  least  one 

alternating  cycle  (Fig  7)  of  odd  length  implying 
that  we  have  at  least  one  PEN,  thus  reducing  the 
problem  to  an  n-1  node  problem. 


PEN  >  ( 


Fig  7 


-7- 


4 .  THREE  ARC  THEOREMS  AJC)  PROOFS 

In  order  to  proof  certain  theorems  about  three  arc 
assignment  solutions  we  introduce   a  third  level  of 
commodities( variables )   kpqrvw 

Aij 
with  similar  types  of  constraints  as  for  the  the  other 
commodities  exept  that  the  quantities  leaving  the  source  and 
transiting  each  node  are  now  equal  to: 

kpqrvw      kpqr 
SUM(j)  Aij      =   Avw. 


kpqrvw      kpqrvw 
3.1  Continuity  at  i:  SUM(j)  {  Aji      -   Aij      1=0 

for  i  other  than  s  or  t; 


kpqrvw   kpqr 

3.2  Flow  from  s  to  k:  SUM(i)  Asi    =Avw, 

kpqrvw    kpqr 
SUM(i)  Aik      =Avw; 

kpqrvw  kpqr 

3.3  Capacity:  0  <=   Aij       <=  Aij, 

kpqrvw  kpvw 

0  <=   Aij       <=  Aij  , 

kpqrvw  qrvw 

0  <=   Aij       <=  Aij , 

kpqrvw  qrkpij 

3.4  Symmetry:  Aij      =  Avw, 

etc.  for  all  permutations, 


Theorem   7: If  we  have  an  RFP  of  length  m<n  then  we 
also  have  an  RFP  of  length  m+1. 

Proof      : Consider  the  graph  of  Fig   9. 

Assume  none  of  the  three  arcs  has  all 

kp 
A   >0  for  (k,p)  on  the  RFP  from  s  to  k. 
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kpqrvw 
This  would  prohibit  commodity  A 
transiting  node  h,  thus  contradicting 
the  flow  constraint  for  this  commodity; 


from 


(s) 


Theorem 
Proof 


kpqrvw     kpqr 
Ah j       = Avw ; 


kp 
Ahd=0 


vw 
Ahe=0 


vw  6.RFP 


Fig  9 


8:  The  RFP's  can  only  be  assignment  solution. 
:The  proof  is  similar  to  the  proof  of  the 
theorem  for  two  arcs . Consider  an  RFP  path  from 
s  to  t  that  does  not  belong  to  an  assignment 
solution.   Due  to  the  flow  constraints  (1.2) 
this  implies  that  we  must  have  at  least  one  arc 
on  one  of  the  remaining  nodes  l.Due  to  theorem 
9  this  node  requires  at  least  one  arc  with  the 
same  RFP  as  on  the  path. Since  it  cannot  go 
back  to  the  RFP  it  has  to  form  an  RF 
circuit, etc . . . 


Theorem   9: A  solution  satisfying  the  AFP  constraints, 
assuming  that  it  is  a  convex  combination  of 
at  most  three  RFP's, must  also  be  a  convex 
combination  of  three  hamiltonian  paths. 

Proof      :It  was  shown  above  that  we  can  only  have 
RFP  'S  along  assignment  solutions.  If  the 
solution  is  not  a  convex  combination  of 
hamiltonian  paths  then  this  would  violate  at 
least  one  of  the  flow  constraints  (1.2). 
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5 .APPLICATIONS 


The  above  formulation  has  been  applied  to  a  limited 
number  of  problems  and  the  initial  computational  results  are 
encouraging. 

I . Peterson  Graph 

Applying  the  above  constraints  to  the  Peterson  graph 
(Fig. 10)  we  can  proove  with  a  relative  small  linear  program 
(  64  variables  ,54  constraints)  that  this  graph  does  not  have 
a  hamiltonian  path  in  its  support,  since  it  does  not  satisfy 
the  AFP  constraints. 


(1)  s,t 


Fig  10 


(1,2) 
Consider  commodity  (1,2)  and  assume  that  A23    >0. 


(1,2)(2,3)    (1,2)(2,3) 
This  requires  that  commodity  A38  +A34       >0. 
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ci,2)(2, 3) 
Assume  A38         =alpha.  This  implies  that   for  each 

(1,2)(2,3)(3,8) 
node  j  we  have      SUM(i)    Aij  =alpha. 

Maximizing  alpha  it  was  found  that  its  optimal  value 

is  zero. 

(1,2)(2,3) 
The  same  holds  for  the  commodity  A34         as  well 
as  for  the  commodities 

(1,2)(2,7) 
A  .Thus  X12  has  to  equal  zero  and  sincethe  graph 

is  symmetric  we  can  conclude  that  we  also  must  have  X15=0. 

Since  a  feasible  solution  to  the  subtour  elimination 
constraints  would  require  that   X12+X15>0  we  have 
a  contradiction  and  thus  proof  that  the  Peterson  graph  does 
not  have  a  hamiltonian  path  in  its  support. 

This  result  is  encouraging  since  most  known  linear 
programming  approaches  to  the  Traveling  Salesman  Problem  have 
failed  at  the  Peterson  graph. 


II .A  six  node  problem (Murtvf  71 ) 

The  traveling  salesman  problem  in  [7] 
following  distance  matrix: 


has  the 


1 
2 
3 
4 
5 
6 


1 

2 

3 

4 

5 

6 

- 

27 

43 

16 

30 

26 

7 

- 

16 

1 

30 

25 

20 

13 

- 

35 

5 

0 

21 

16 

25 

- 

18 

18 

12 

46 

27 

48 

- 

5 

23 

5 

5 

9 

5 

- 

As  reported  in  [3]  the  optimal  assignment  solution  has  a 
value  of  54  units. The  optimal  solution  using  the  subtour 
elimination  constraints  (Fig  11)  has  a  value  of  57.9.   Using 
a  branch  and  bound  procedure  the  optimal  tour  can  be  found 
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Fig. 11 


From  Fig.  11  we  notice  for  example  that  it  is  not  feasible  to 
send  0.7  units  of  commodity  X35  from  s  to  node  3,  and  0.7 
units  of  slack  from  node  5  to  t . 


Introducing  flow  requirements  for  commodity  A35: 

(3,5) 
i.e.  SUM   Xsi     =X35, 
i 

(3,5) 
SUM   Sit      =X35, 
i 

Continuity  at  node  i : 


35  35  35 

SUM     Sij  -   SUM   Sji=0,    SUM   Xij 
J  J  3 


35 
SUM   Xji  =0, 
J 


35  35 

SUM   Sij   +  SUM    Xij   =X35  ; 
J  J 

The  linear  programming   solution  is  integer  valued  and  is 
thus  the  optimal  tour ( Fig. 12 ). It  has  a  length  of  63  units. 
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s,t  (1)^ 


Fig. 12 


6 .PROPOSED  ALGORITHM 


As  mentioned  above,  numerical  results  are  quite  limited  at 
this  time,  but  the  following  algorithm  appears  to  be 
promising: 


STEP  1   Find  the  optimal  solution  using  the 

assignment  procedure. If  we  have  a  tour  STOP, if 
not  go  to  STEP  2. 

STEP  2   Find  an  optimal  solution  using  the  two-commodity 
subtour  elimination  constraints  [3]. 
If  integer  valued  STOP, if  not  go  to  STEP  3. 

STEP  3   Find  an  arc  with  violated  flow  requirments; 

(If  none  can  be  found  go  to  STEP  4)  introduce  the 
necessary  constraints  and  find  the  optimal 
solution. If  the  solution  is  integer  valued  STOP 
the  optimal  tour  has  been  found,  else  we  have 
a  lower  bound. Introduce  this  lower  bound  (i.e. 
the  objective  function  value  has  to  be  >=  this 
bound)  and  discard  the  previous  flow 
requirements.   Go  to  STEP  3. 

STEP  4   Find  a  two  -arc  commodity  with  violated  flow 
requirements 

STEP  k-i-2  Find  a  k  arc  commodity  
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7 .CONCLUSIONS 

Although  the  number  of  potential  variables   is  extremely 
large  the  above  outlined  algorithm  would  only  need 
4-commodities  at  a  time  (  and  is  thus  computationally 
feasible  even  for  very  large  problems),  if  at  each  STEP  we 
get  an  increase  in  the  objective  function  value  (No 
alternative  optimal  solutions)  and  probably  not  too  many 
commodities  if  we  do  have  alternative  optimal  solutions. 


At  each  STEP  we  introduce  subtour  elimination  constraints 
for  the  commodities  introduced  in  the  previous  STEP. 
The  interesting  theoretical  question  arises  as  to  wether  we 
have  a  reasonable  worst  case  upper  bound  on  the  number  of 
STEPS  we  do  have  to  execute. If  this  upper  bound  would  be  a 
constant  then  this  would  imply  that  P=NP.Thus  it  appears 
likely  that  this  upper  bound  is  an  exponential  function  of 
the  number  of  nodes  and  or  arcs  of  the  graph. 


Past  computational  experience  with  the  two-coirunodity 
approach  [3], [5]  has  shown  that  the  larger  the  problem  size 
the  more  likely  the  optimal  solution  to  the  two- 
commodity  problem  will  be  integer  valued;  thus  we  may 
expect  that  for  large  problems  we  rarely  need  more  than 
4-commodities  and  that  in  view  of  the  results  for  the 
Peterson  graph  we  probably  dont  have  to  go  very  often  beyond 
STEP  4  or  5  . 

Extensive  numerical  experimantation  is  under  way   at  this 
time,  and  will   be  reported  on  soon. 
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